Fractal design for an efficient shell strut under gentle compressive loading 
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Because of Euler buckling, a simple strut of length L and Young modulus Y requires a volume of 
material proportional to L 3 / 1 / 2 in order to support a compressive force F, where / = F/YL 2 and 
/ < 1. By taking into account both Euler and local buckling, we provide a hierarchical design for 
such a strut consisting of intersecting curved shells, which requires a volume of material proportional 
to the much smaller quantity L? f exp \2^J (In 3)(ln • 

PACS numbers: 46.25.-y, 46.70.De 



I. INTRODUCTION 

Fractals [lj occur ubiquitously in nature and appear to 
arise in many different ways; examples from the physical 
and biological sciences include colloidal nocculation@, Q , 
percolation phenomenal and the structure of transport 
networks in organisms [5j. In the area of structural me- 
chanics, it has been claimed that the fractal morphology 
of trabecular bone is responsible in part for its mechanical 
efficiency @; while the complex suture patterns of am- 
monites of the Jurassic and Cretaceous periods have been 
conjectured to give greater strength to their shells@, Q. 

Recent theoretical work on a highly simplified model 
system consisting of a brittle pressure-bearing plate has 
also suggested that a fractal structure can be highly ef- 
ficient when the loading conditions are very gentle and 
the material very brittle [§] . It is therefore of interest to 
explore other circumstances under which fractal design 
principles can lead to high mechanical efficiency, in case 
there is a general theorem underlying the optimal design 
of elastic structures under gentle compressive loading. To 
this end, we consider here the buckling behaviour of com- 
pression members. 

Consider first the classical case of an Euler strut [T3| 
in the form of a solid, cylindrical column of radius ro and 
length L, made from an isotropic, linear elastic material 
of Young modulus Y and Poisson ratio v, and subject to 
a compressive force F at the freely hinged ends. 

We define two non-dimensional parameters: / = 
F/(YL 2 ), which is the compressive force scaled by YL 2 , 
and v = tttq/L 2 , which is the volume of material used, 
scaled by L 3 . We are interested in the regime of gentle 
loading, by which we mean / < 1 and »<1. 

Because of Euler buckling, the strut can only with- 
stand forces such that F < n 2 YI /L 2 , where Iq = ^r^/A 
is the second moment of the cross sectional area about 
the neutral axis of the beam [l(| [H| • Therefore the (non- 
dimensionalised) volume of material required to with- 
stand a load represented by / is given by 



»(0) = 27T- 1 / 2 / 1 



/2 



(1) 



where we have neglected the material required to make 
the freely hinged couplings at the ends of the strut. For 
the purposes of this paper, we call the above solid strut 
a "generation G — 0" structure, and this is written as an 
argument for the volume variable v in Eq. (TTJ). 

We note in passing that the contrast in efficiency be- 
tween compression members and tension members (for 
which v oc /) is a persistent theme in structural engi- 
neering. The scaling of Eq. ((T|) plus the cost of cou- 
plings mean that efficient structures tend to have few 
compression members and many long tension members; 
the paradigmatic example being a tent [H, EH • 

To define a, G — 1 structure, we choose a hollow cylin- 
drical shell .which is also a classic problem in elasticity 
theory [HI EH- We choose the length as always to be L, 
and we denote the radius by r% i where the first index 
refers to the "generation number" G = 1 and the second 
index will be explained when we describe structures of 
higher generation number. This G = 1 structure consists 
of one cylinder parallel to the compression direction, and 
we express this trivial fact by = 1. The thickness of 
the sheet of elastic material making up the curved surface 
of the cylindrical shell is denoted by si t i, which specifi- 
cally represents the volume of material required to make 
unit area of the curved surface. We call this quantity the 
"material thickness" of the curved sheet. We also intro- 
duce an "effective elastic thickness" ti i for the curved 
sheet. For the generation 1 structure, the curved sheet 
is simple in topology and uniform in thickness, and so 
si i = ii i. Lastly we have an effective Young modulus 
Yi.i and Poisson ratio v\ t \ for the sheet. For a G = 1 
structure Yj. i = Y and i^i.i = v. In all of these expres- 
sions, the first index refers to the generation number, and 
the second index will take values from 1 up to the gen- 
eration number of the structure, as will be explained in 
section IIIIl 

Provided that t\^\ <C fi i <C L, the column now has a 
second moment of cross sectional area about the neutral 
axis given by I\ = irsi^rf 1 and the volume of material 
used to construct it is given by v — 27rri i iSi i i/L 2 . The 
requirement that Euler buckling not occur then imposes 
the constraint F < ir 2 Yx t iIi/L 2 or 



'Electronic address: robert.farr@unilever.com 



f>2/ 1 / 3 (t 1A /L) 



2/3 



(2) 
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In contrast to the solid column, there is now the pos- 
sibility of local buckling. This happens when [TJ llij 



F 



27^1,1%! 

3(1 



(3) 



which provides a second constraint on the structure re- 
quired to support a force F. 

The generation 1 structure with the highest mechanical 
efficiency (smallest value of v for a given /) is therefore 
specified by: 



«(1)=2 



3(1 - v 2 



4ir 2 



1/6 



f2/3 



with 



3(1 - 



4tt 2 



-1/12 



f l/6 



«1£ 



3(1 -v 2 ) 



Air 2 



1/4 



1/2 



(4) 



(5) 



(6) 



and therefore t\,x -C 7%! <I as assumed above. 

Eq. (j4|) represents a considerable gain in efficiency over 
the solid column of Eq. ([T]), but does not rival the effi- 
ciency (v oc /) which can be obtained for tensile loading. 

We note that both here and in all subsequent sections, 
we have taken a conservative approach in calculating v, 
by assuming that the structure fails when the first buck- 
ling bifurcation is encountered. Engineering structures 
(especially shells) will often support considerably higher 
loads in the post-buckling regime before catastrophic fail- 
ure [ll|. Although such complexities are certainly of 
practical importance, we choose to ignore them in this 
investigation and try, where possible, to obtain estimates 
which are upper bounds for v. 



II. A COMPOSITE PLATE 

When designing a plate or shell which may buckle, it 
is standard engineering practice to introduce stiffening 
plates, longitudinal stringers, bulkheads, or similar de- 
vices in order to stiffen the structure and/or suppress 
buckling modes [U, [l3| . 

In this paper we take a similar approach, but re-design 
the structure in a systematic and hierarchical manner, 
which can be iterated in the limit / — > 0. We do 
not presume to do this in the most efficient manner 
(we are almost certainly over-engineering the protection 
against many of the buckling modes we wish to avoid) 
but nevertheless the design we describe allows us to sys- 
tematically change the scaling of v with / and there- 
fore to approach more closely the scaling which can be 
achieved for a rod under tension; achieving ultimately 



v oc / exp 



2 v /(ln3)(ln/- 



In the limit / — !• 0, this is 



smaller than any scaling of the form v oc f" with f3 < 1. 

To proceed in this direction, consider first of all a sim- 
ple thin, fiat plate of uniform thickness t, lying in the x—y 
plane and made out of an isotropic elastic material of 
Young modulus Y and Poisson ratio v. Suppose further- 
more that this plate may be deformed by applied stresses, 
leading to stretching or shearing of the middle plane (l5j. 
and also possibly to out-of-plane deflections which may 
be large compared to the plate thickness. An appropriate 
approximation to describe the behaviour of the plate is 
due to von Karman 16]. In this theory, the equilibrium 
behaviour (given suitable boundary conditions) can be 
obtained by minimising an energy functional. The elas- 
tic part of the energy (as opposed to that from external 
forces) can be written as the sum of two terms [IT}: the 
energy Us associated with stretching of the middle plane 
of the plate, and a bending energy Ub- 

If the 2-dimensional strain tensor for the middle plane 
of the plate is given by e[x, y), then the stretching energy 
will be given by [3, [TJ [lT} 



Us 



Yt 



2(1 - £> 2 



dxdy jz> [Tr(e)] 2 + (1 - £>)Tr(e 2 )| . 

Furthermore, if the plate is bent out-of-plane in the z- 
direction, by an amount wix, y), then the bending energy 
stored will be given by [17[ 

Ub = 2A ^~ 2) J dxdy {[Tr(iT)] 2 - 2(1 - D)det(H)} , 

(8) 



where H is the Hessian matrix 

/ d 2 w d 2 w 



H{x,y) 



dx 2 dxdy 
d 2 w d 2 w 



(9) 



dxdy dy 2 



Now consider the composite plate shown in the bot- 
tom part of Fig. [T] This structure is built from three 
intersecting sub-structures (as shown separately in the 
top three parts of Fig. [IJ , each identical, save for being 
rotated ±27r/3 radians relative to one another. Where 
the sub-structures pass through one another, we imagine 
them being joined or welded along their curves of inter- 
section. Each sub-structure consists of an infinite set of 
parallel hollow right circular cylinders, with their axes 
in the x — y plane, and placed so that each touches two 
neighbours and is welded to each of its two neighbours 
along these lines of contact (which also lie in the x — y 
plane). We specify all these welds so we can be sure that 
on deformation at long length scales, the plate behaves as 
a single entity, rather than separating into its constituent 
cylinders. 

Each of the component cylinders has a radius r and a 
wall thickness t <£L r. Because we have chosen the com- 
posite plate to have six-fold rotational symmetry about 
the z-axis, then on long enough length scales the compos- 
ite plate must behave elastically as though it is isotropic 



3 





FIG. 1: Bottom image is part of a composite plate, which 
is constructed by merging the three substructures shown in 
the top three images of the figure. Each substructure is an 
infinite set of parallel right circular cylinders which are joined 
along their lines of osculation. The three sub-structures are 
identical save for being rotated by ±27r/3 radians relative to 
one another about an axis perpendicular to the plane. Each 
right circular cylinder has a radius r, a wall thickness t <C r 
and is infinitely long. 



under rotations about the z-axis. This is because to lead- 
ing order, the elastic stiffness of the plate under stretch- 
ing and bending is represented by rank 4 tensors in two 
dimensions (which are contracted with two dimensional 
rank 2 deformation tensors to form the scalar energy). 
These rank 4 elastic tensors may be invariant under ro- 
tation about the z-axis (and so be consistent with any 
rotation group C„), or may have one of the symmetry 
groups Ci or C4. However the composite plate we have 
described has the symmetry group Cq, which is not con- 
sistent with C2 or C4. 

Calculating the stretching and bending energies of the 
composite plate shown in the bottom panel of Fig. [T] is 
non-trivial even for the long-wavelength (much greater 
than r) deformations in which we are interested. In- 
deed, obtaining the correct numerical pre-factors would 
require a finite element calculation for the structure. In- 
stead, and in order to proceed rigorously and in an ana- 
lytical manner, we make an approximation which under- 
estimates the stiffness (over-estimates the compliance) of 



the structure, so that we would end up eventually with 
an upper bound on the minimum value of v required for 
the final columns in the sections to follow. 

The approximation we make is a "ghost approxima- 
tion" , in which we imagine that the cylinders composing 
the composite plate of Fig. [T] are no longer welded to- 
gether, but are free to move past and indeed through one 
another; but all follow the imposed deformation field. 

Under this approximation, consider what happens to 
one of the constituent cylinders when the composite plate 
is subjected to an in-plane stretching deformation, rep- 
resented by a two-dimensional strain tensor e(x, y) with 
principal components e\ and e%. Let the cylinder be at 
an angle relative to the direction of the principal com- 
ponent e\. The cylinder will then be stretched parallel 
to its symmetry axis with a strain 



r](6) = e\ cos 2 9 + e2 sin 2 9. 



(10) 



It may also rotate, but because of the ghost approxima- 
tion, it experiences no resistance to this motion. 

Adding up the stretching energies for the three sub- 
structures, we find a lower bound for the total stretching 
energy of the composite plate Us.com given by 

£/s,com > ^ f dxdy {3 [Tr(e)] 2 + 6Tr(e 2 )} . (11) 

Consider next a single cylinder when the composite 
plate is subjected to an out-of-plane bending deforma- 
tion field w{x,y), which is slowly varying in x and y 
(compared to the cylinder radius). If the cylinder is at 
an angle 9 to the x-axis, then it will have an out-of-plane 
curvature given by 



k{8) = 



+ 



d 2 w . „„ d 2 w . 9 
■ sin 29 + -rr^r Sin ' 



(12) 



dx 2 dxdy dy 2 

From thin-beam theory [ll| it will therefore have an elas- 
tic energy per unit length given by 

1, 



u{9) = -Yixtr 6 [k(6)Y 



(13) 



Adding up contributions from the three sub-structures, 
the bending energy C/B iCom of the entire composite plate 
therefore has a lower bound given by 



U B , 



> 



3irYr 2 t 



32 



dxdy {3 [Tr(iJ)] 2 - 4det(#)j 



( 14 ) 

If we compare the lower bounds on elastic energy from 
Eqs. (fTTj) and (|T4|) with those for a uniform plate (Eqs. ([7]) 
and ) then we see that for long wavelength defor- 
mations, the plate is at least as stiff as a uniform plate 
with effective thickness, Young modulus and Poisson ra- 
tio given by 



Y e 



cir 



i/ e s = 1/3. 



(15) 
(16) 
(17) 
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Furthermore, the composite plate uses an amount of 
material per unit area given by 

s cff = 37rf. (18) 

The results in Eqs. (|T5|) , (H~6l) and (|TT|) provide us with 
the information required to calculate lower limits on the 
stresses required to produce buckling on length scales 
much larger than r. However, a composite plate of this 
kind can also fail through local buckling by one or more 
of the constituent cylinders undergoing a local buckling 
instability. 

This can be dealt with analytically for an isolated 
cylinder, which is the result used above for local buckling 
of a generation 1 structure (Eq. (j3J). However, without 
an extensive finite element study, covering a range of pa- 
rameters, it is much harder to provide a good lower bound 
on the in-plane compressional stress required to excite 
these modes in the composite plate. This is because the 
stresses in intersecting cylinders could potentially gener- 
ate buckling, rather than having no effect (as in the ghost 
approximation) or suppressing these modes. 

In what follows, we will make the simple, but not nec- 
essarily accurate approximation that provided the largest 
compressive principal stress component lies in the direc- 
tion of the axes of one of the substructures in the com- 
posite plate, then we will get local buckling only under 
the same circumstances as for an isolated cylinder. This 
is also a kind of "ghost" approximation, in that the role 
of the other parts of the substructure are ignored, but 
although it is plausibly a conservative approximation, it 
no longer provides a strict bound on v. 

III. THE GENERATION 2 STRUCTURE 

Imagine taking the curved cylindrical shell which forms 
the generation 1 structure of section Q] and replacing the 
solid curved shell with a composite plate similar to that 
in Fig.[TJ but curved to follow the original cylindrical sur- 
face. An example is shown in the bottom image of Fig. [21 
with the top three images in the same figure showing the 
three substructures which are merged to form the final 
column; one of the substructures has the axes of its con- 
stituent cylinders aligned with the long axis of the col- 
umn itself. We refer to the result as a generation G = 2 
structure. 

The length of the generation 2 structure is taken to 
be L, and the radius of the entire column is f2,2- At 
the largest scale, there is only ever one cylinder, and we 
represent this trivial fact by n.2,2 = 1- However, that 
substructure making up the composite shell which has 
cylinders aligned with the column length is composed of 
more than one cylinder, and we denote this number by 
712,1. For the structure in Fig. [21 we have ri2,i = 12, 
as can be seen more clearly in the section through the 
relevant sub-structure shown in Fig. [31 

The thickness of the thinnest shells, which make up the 
cylinders of the composite shell is t2,i, and these have 




FIG. 2: The bottom image is an example of a generation 
2 structure or column, which consists of three intersecting 
sub-structures that are shown in the top three images. Each 
sub-structure is composed of hollow cylindrical shells, with 
one sub-structure having the cylinders aligned parallel to the 
long axis of the entire column, and in the other two sub- 
structures the cylinders are wrapped in a helical arrangement 
(left handed for one sub-structure and right-handed for the 
other) . 



Young modulus Y 2 .\ = Y and Poisson ratio 1/2,1 = v. 
These shells form the cylinders of the composite sub- 
structures, which each have radius r2,i 3> £2,1- The re- 
sulting composite shell has an effective clastic thickness 
f2,2j an d effective Young modulus and Poisson ratio given 
by Y 2j2 and v 2 ,2- 

The geometrical terms are illustrated in Fig. [3j which 
for clarity shows only a cross-section through the sub- 
structure which has its cylinders aligned with the axis 
of the entire composite column. We note finally that in 
the case T2,i *C ^2,2 then we can count the number of 
cylinders in this sub-structure using 

"2,1 =7IT2,2/r2,l. (19) 



Rather than analysing the efficiency of the G = 2 struc- 
ture here, we proceed directly to the general case in the 
following section. 
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FIG. 3: Section through one of the sub-structures in Fig. [2] 
which has the component hollow cylinders parallel to the axis 
of the entire generation 2 structure. The figure shows defi- 
nitions of various lengths required to specify a generation 2 
structure. 



IV. THE GENERATION G STRUCTURE 

To make the generation 3 structure, we imagine re- 
placing all the curved but solid shells in a generation 2 
structure (which are of uniform thickness £2,1) with com- 
posite shells, as described for a flat plate in section [TT1 
above. 

The thinnest (solid) shells comprising this new struc- 
ture have a thickness of ts,i and compose thin cylinders of 
radius r^ t \. These form the substructures of curved shells 
with effective elastic thickness £3,2- These curved shells 
form cylinders of radius which compose the substruc- 
tures of a new composite shell, which is of effective elastic 
thickness ^3, and from this final doubly-composite shell, 
a hollow cylinder of length L and radius 7-3,3 is formed, 
which is the final generation 3 structure. 

By iterating this process, we end up with a genera- 
tion G structure. In general, r^m is the radius of each 
of the (usually composite) cylinders at hierarchical level 
m of the generation G structure: Tg,i is the radius of 
the smallest (non-composite) cylinders and tq,g that of 
the multiply composite column itself. Similarly ^G,m is 
the effective elastic thickness of the curved (and usually 
composite) shells making up the cylinders at hierarchical 
level to in the generation G structure. 

We assume that 

tc,i < r G ,i < t G>2 < r G , 2 < . . . < t G ,G < f G ,G < L 

(20) 

and by definition we can say the following: 



*G,1 = Sg,1, 

Y G ,x = Y, 

VGA = V. 



(21) 
(22) 
(23) 



are related to one another through Eq. (fT8|) by 

S G ,m = 37TSG,m-l 



(24) 



and therefore the total (non-dimensionalised) volume of 
material used is given by 



v(G) = 2nLrG.G s G,G / 



(25) 



The effective clastic properties for in <E {2, . . . , G} are 
related to one another in the ghost approximation (which 
gives an upper bound on v) by results analogous to Eqs. 
(1T51) . (|To| and <[TTJ) of section QTJ namely 



Yt 



G,m 



= ^ ( ^=L\ Y G<m -i, 

m = 1/3. 



VQ \rG,m-l 
VG 



(26) 
(27) 
(28) 



At the largest length scale (i.e. the column itself), 
the structure is subject to Euler buckling and therefore 
the largest (non-dimensionalised) force it can support is 
subject to the constraint 



/< 



IG*G,G 



YL 2 



(29) 



where the relevant second moment of the cross sectional 
area about the neutral axis is given by 



I G = xtG,Gr G ,G- 



(30) 



For the constraints due to localised buckling, we pro- 
ceed as follows: at the largest length scale, there is one 
(multiply composite) cylinder aligned with the long axis 
of the column. This fact is captured by the equation 
nG,G = 1. 

At the other length scales, we count the number of 
smaller cylinders in the substructures which are aligned 
with the long axis of the entire column in the following 
recursive manner, based on Eq. (|19|) : 



(31) 



rG,m 
ro.rn-1 



where m g {2, . . . , G}. 

Each of these cylinders at a level m of the structure has 
a shell with an effective elastic thickness of tG, m , effective 
Young modulus YQ,mi effective Poisson ratio vc.m and 
supports a force no more than 



F G,r 



(32) 



The material thickness of the different composite shells 



which we obtain by neglecting the support provided by 
the other two helically arranged sub-structures at this 
level (and so on recursively). 

As discussed in section [IT] we make the crude approxi- 
mation that the local buckling condition for a composite 
shell can be obtained from that for the isolated cylinders 
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composing it. This again uses the ghost approximation, 
but in this case the approximation no longer provides 
a strict bound. The result is a sequence of conditions 
to avoid local buckling at each hierarchical level in the 
structure, analogous to Eq. ([3]) and given by 



2TTY G ,mt 2 G , 



3(1 J 



(33) 



We now proceed to solve the recursion relations of 
Eqs. flUHM}, flUH} and d3TJ) , keeping first of all t GA 
and rQ^ m as parameters for optimisation: 



t G 



Y G , m = 



SG,rn = (37T) m t G ,\ 

tc,i m = 1 
V6ra, m -i 2 < m < G 

Y TTl — 1 

1I ^_(_ta^_\ Y 2 <m<G 

{v m = 1 
1/3 2<m<G 

G-m r G,G 



n G , 



(34) 
(35) 

(36) 

(37) 
(38) 



The condition to be at the limit of Eulcr buckling 
(Eq. ([29])) therefore becomes 



f 



r G+2 



' G.G 



tGA, 



(39) 



the condition to be at the limit of local buckling at the 
smallest scale of the structure is (from Eqs. (I3"2"j) . (|3"3")l 
and (38])) 



/ 



2tt g 



r G,Gt G ,l 



(40) 



and to be at the limit of local buckling at the other levels 
in the structure, gives (from Eqs. ([52]). (|3"5"]) . and 
(US}) for 2 < m < G 



/ = 



37T G TQ m —\rQ QtQ i 



L- 



r G, 



(41) 



Eqs. (HI]), (ED), (32]), & and (H]) can be solved for 
G > 1 to give finally 

v(G) = 2( 1 + G )/( 2 + G )3( 2G2 - 1 )/[ 2 ( 2 + G )] 7 r( G - 2 )/( 2+G ) 

X(l - i y 2 ) 1 /[2(2+G)] jl-[l/(G+2)] ; ^ 

t"G,l = 2 ^3/[2(2+G)] 3 (9-6G)/[4(2+G)] 7r -(2G 2 -G+2)/[2(2+G)] 

X(l - ^2)3/14(2+0] j(l+2G)/[2(2+G)^ 43 ) 

r _GjG = 2 l/[2(2+G)] 3 (2G-3)/[4(2+G)] 7r -(3G+2)/[2(2+G)] 

L 

X(l - ^ 2 )-V[4(2+G)] yl/[2(2+G)] j( - 44 > ) 



TABLE I: Example calculation for the mass M if a structure 
required to support F = lOkN over a distance of L — 200m 
when the structure is made from a material similar to steel, 
with Y = 210GPa, v = 0.29 and density p = 8000kgm -3 . 



G 


M 


ta,i 


ra,i 


rc,2 


ra,3 





79 tonnes 




12.5cm 






1 


941kg 


0.12mm 


81cm 






2 


319kg 


1.4/xm 


1.2mm 


2 Am 




3 


261kg 


63nm 


14/xm 


8.2mm 


4.6m 



O 




-20 -18 -16 -14 -12 -10 

log 10 (/) 



FIG. 4: Plot of log 10 (u) versus log 10 (/) for the optimal 
structures with G = 0,1,2,3,4,5,6, using Eq. (|42|) from 
the text and v = 0.29. The curves are only drawn for the 
range of / where rc,a/L < 1/20, ra, m /rG,m+i < 1/20 and 
tG,i/ra,i < 1/20, which we take as an approximation to the 
condition of Eq. (1201) in the text. 



while for 1 < m < G 



2(2G-2m+l)/[2(2+G)]g-(12G-14m+3)/[4(2+G)] 
-(2G 2 -2Grn-G+4m+2)/[2(2+G)] 



X7T 



x (l _ , y 2)-(2G-2m+l)/[4(2+G)] yr(l-2fn+2G)/[2(2+G)] Jax\ 



As a simple practical example, consider a strut of 
length L = 200m which is required to support a force 
of F = lOkN, and which is made from a model mate- 
rial, similar to steel, with Y = 210GPa, v = 0.29 and a 
density of 8000kguT 3 , so that / = 1.2 X 10~ 12 . 

A cable supporting this force under tension would re- 
quire a mass of 8kg (assuming a yield stress for the mate- 
rial of 200MPa, and neglecting the mass of couplings at 
the ends). The masses (M) of "steel" required for various 
structures described in this paper are shown in Table Q] 

Lastly, we note that for a given value of /, several gen- 
erations of structures may be compatible with the con- 
ditions of Eq. (|20l) , as illustrated in Fig. @] In the limit 
/ —> 0, we can calculate the envelope of these curves in 
order to obtain the global optimally efficient structure 
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within this class, through solving [R 

v = v(G) 
dv(G) 



dG 



0, 



(46) 
(47) 



where v(G) is given by Eq. (|42j) . 

In the limit of small /, we can expand the exponent 
of Eq. (|4"6")) in powers of 1/G to obtain the assymptotic 
results 



2tt/ 



exp 



2V(ln3)(ln/-i 



ln3 



(48) 
(49) 



CONCLUSIONS 



We have described compression members consisting of 
intersecting curved shells in a fractal or hierarchical ar- 
rangement which are highly mechanically efficient in the 
limit of light compressional loading. 

Fractal designs for efficient plates under gentle pres- 
sure loading have recently been studied in Ref. [9j. In 
this work, the fractal design arises from two competing 
tendencies in the structure: Firstly there is a geometri- 
cal feature of the plate (narrow, tall spars) which when 
developed to extremes can lead to very high mechani- 
cal efficiency. Secondly, there is a limit to how far this 
feature may be developed, which is ultimately a vulner- 
ability to buckling. One spar can however be used to 



provide partial support for another, and this leads to the 
final hierarchical design. 

The parallels with the problem of the present paper 
should be apparent, and so we suspect that fractal forms 
may be a general property of optimal elastic structures 
under gentle and at least partially compressive loading. 

At this stage, we are not able to frame a precise math- 
ematical conjecture, and we do not know how other pa- 
rameters will figure in the analysis. In previous work brit- 
tleness Q was important, and it seems highly probable 
that in the present work, the amplitude of imperfections 
in either the geometry or the uniformity of the loading 
could be crucial to determining the mechanical efficiency 

We therefore hope that further examples, and perhaps 
ultimately theorems, will shed light on a problem whose 
geometrical solutions promise to be useful and even beau- 
tiful. 
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